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WlO .  Abstract. 
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the  free  stream  streamsheet  due  to  the  blade-to-blade  velocity  gradient,  the  magni  - 
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perature  rise  are  similar  to  experimentally  measured  values.  Thus,  this  calculation 
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ABSTRACT 


Current  through-flou  calculation  schemes 
provide  the  compressor  designer  with  a  2-dimensional 
inviscid  flow  analysis  of  a  blade  row.  Due  to  blade 
curvature  and  viscous  effects  the  flow  field  is  quite 
3-dimensional.  A  numerical  scheme  has  been  developed 
to  calculate  the  radial  migration  of  secondary  flows 
in  a  rotating  blade  row  and  the  influence  of  these  flows 
on  the  energy  transfer  of  the  rotor.  A  simple  set  of 
equations  far  the  radial  equilibrium  and  resulting  tra¬ 
jectory  of  these  fluid  particles  is  derived.  A  boundary 
layer  calculation  was  added  to  correctly  model  rotor 
boundary  layer  fluid.  The  calculation  scheme  was  then 
applied  to  the  2-dimensional  flow  field  obtained  from 
an  existing  through-flou  analysis  of  a  transonic  rotor. 

The  results  showed  the  twisting  of  the  freestream  stream- 
sheet  due  to  the  blade- to-blade  velocity  gradient,  the 
magnitude  of  the  boundary  layer  migration  and  the  in¬ 
fluence  of  the  radial  shift  on  the  enthalpy  rise  of  these 
particles.  The  resulting  values  of  radial  velocity  and 
temperature  rise  are  similar  to  experimentally  measured 
values.  Thus  this  calculation  scheme  can  provide  the  com¬ 
pressor  designer  with  a  quick  first  approximation  of  the 
influence  of  3-dimensional  and  viscous  effects  on  the 
flow  field  of  a  rotating  blade  row. 
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1 .  INTRODUCTION 

Current  compressor  analysis  tools  include 
through-flou  calculation  schemes  which  provide  the  com¬ 
pressor  designer  with  an  inviscid  axisymmetric  view  of 
the  flow  field  within  a  blade  row.  Such  results  are 
useful  for  determining  the  radial  variations  in  the  flow 
field  for  a  mean  blade- to-blade  streamline.  Such  a 
meridional  streamline  network  is  shown  in  Figure  1  for 
a  transonic  rotor.  Blaae  curvatures  and  viscous  effects 
however  cause  considerable  distortion  of  this  flow  field. 
The  blade  curvatures  produce  a  velocity  graoient  in  the 
blade- to-blade  plane  creating  a  non-axisymmetric  flow 
field.  This  velocity  gradient  will  also  cause  some 
warpping  or  twisting  of  the  blade- to-blade  streamsheet 
(  see  Figure  2.  ).  Viscosity  of  the  fluia  will  produce 
hub  and  tip  casing  and  blade  boundary  layers  and  wakes 
producing  a  highly  three-dimensional  flow.  Thus  while 
being  a  useful  engineering  tool  the  through-flow  calcu¬ 
lation  needs  to  be  supplimented  with  a  procedure  for 
determining  these  additional  three-dimensional  aspects 
of  the  flow  field. 

Due  to  the  fluid  viscosity  some  of  the  fluid 
particles  in  the  flow  field  will  have  a  deficit  or  defect 
in  velocity  relative  to  the  essentially  inviscid  free- 
stream  flow.  These  low  momentum  fluid  particles  are 
subjected  to  the  same  pressure  forces,  both  circumferen- 
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tially  and  radially,  as  the  freestream  but  due  to  their 
lower  momentum  their  path  or  trajectory  through  the  blade 
row  will  be  influenced  more  by  these  forces.  This  momen¬ 
tum  deficit  is  thus  the  source  of  some  of  the  secondary 
flows  in  the  blade  row.  A  classical  example  of  such  sec¬ 
ondary  flows  is  the  migration  of  endwall  boundary  layer 
fluid  from  the  pressure  surface  to  the  suction  surface 
of  a  turbine  nozzle  uane  row.  The  radial  migration  of 
rotor  blade  boundary  layer  fluid  is  another  example  of 
these  flows. 

Several  experimentam  investigations  have  pre¬ 
sented  complete  three-dimensional  measurements  of  a  rotor 
exit  flow  field  (REF  1,2, 3, 4).  These  studies  have  shown 
that  radial  velocities  in  a  rotor  wake  (due  to  the  radial 
migration  of  boundary  layer  fluid)  can  be  quite  large. 
Values  of  these  measured  radial  velocities  range  from 
10^  -  1 5%  of  the  axial  velocity  in  a  low  speed  rotor  (REF 
2,3)  to  reportedly  as  much  as  100/b  of  the  axial  velocity 
for  a  transonic  rotor  (REF  4).  It  has  also  been  shown 
that  the  magnitude  of  the  radial  velocity  in  the  wake  is 
influenced  by  the  blade  loading  (REF  3). 

One  aspect  of  these  secondary  flows  that  has 
not  recieved  much  attention  is  how  the  energy  transfer 
of  the  rotating  blade  row  is  influenced  by  these  flows. 

It  was  originally  believed  that  due  to  the  relative 
motion  of  the  rotating  and  stationary  blade  rows  circum- 
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ferential  variations  in  the  flow  field  would  be  quickly 
"mixed  out".  Measurements  (REF  5)  have  however  shown 
that  variations  in  the  total  temperature  as  large  as  10^ 
of  the  overall  stage  temperature  rise  can  still  exist  at 
the  exit  of  a  stator  row  downstream  of  the  rotor.  These 
variations  in  total  temperature  could  only  have  been  caused 
by  the  uneven  energy  addition  of  the  rotor.  It  is  believed 
that  the  secondary  flows  in  the  rotating  blade  row  are 
the  source  of  these  temperature  variations.  The  magni¬ 
tude  of  these  variations  is  significant  when  attempting 
to  calculate  rotor  or  stage  efficiency  from  measured  exit 
total  temperature.  The  existing  temperature  gradient  can 
also  be  the  source  of  thermally  induced  secondary  flows 
(REF  6) . 


Smith  (REF  7)  showed  that  the  energy  transfer 
to  a  fluid  particle  in  a  turbomachine  is  proportional  to 
the  time  the  particle  remains  within  the  rotating  blade 
row.  Thus  low  momentum  particle  that,  due  to  their  velo¬ 
city  defect,  remain  within  the  rotating  frame  longer  will 
acquire  a  larger  enthalpy  increase  than  freestream  fluid. 
These  low  momentum  particles  will  also  have  a  tendency 
to  move  radially  outward  due  to  an  imbalance  of  pressure 
forces  and  centrifugal  forces.  This  change  in  radial  lo¬ 
cation  will  contribute  further  to  the  increased  energy 
transfer  to  the  low  momentum  particles. 


With  a  simple  wake  model  Kerrebrock  and 
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Wikola jczak  (REF  8)  shoued  that  blade  boundary  layer 
fluid  due  to  a  slip  velocity  of  the  resulting  uake  uill 
tend  to  be  separated  from  the  freestream  fluid  as  the 
flow  progresses  through  the  successive  blade  rows.  Thus 
the  low  momentum  but  high  total  temperature  uake  fluia  of 
the  rotor  will  tend  to  be  collected  by  the  downstream 
stator  and  appear  in  the  wake  of  the  stator  blade.  Thus 
areas  of  concentrated  high  total  temperature  rotor  wake 
fluid  can  occur. 

Experimental  evidence  of  circumferential  total 
temperature  variations  in  the  flow  field  has  been  pre¬ 
sented  by  Keenan  (REF  5),  Leboeuf  (REF  9)  and  Kool  (REF 
10).  Analytical  attempts  at  predicting  the  influence 
of  secondary  flows  on  the  energy  transfer  and  the  result¬ 
ing  temperature  field  are  scarce.  Before  such  an  analysis 
can  be  made  it  is  important  to  properly  predict  the  sec¬ 
ondary  flows  themselves.  Thus  a  significant  portion  of 
this  project  is  devoted  to  establishing  a  simple  yet 
physically  correct  model  for  the  low  momentum  particles 
in  a  rotating  blade  row.  Such  a  model  should  include  not 
only  the  influence  of  viscosity  but  also  the  effects  of 
blade  and  meridional  flow  path  curvatures. 

Throughout  this  analysis  the  author  took  a 
"particle"  approach  for  formulating  and  mentally  visual¬ 
izing  the  resulting  secondary  flows.  This  involved  as¬ 
suming  an  inviscid  fluid  particle  that  for  some  reason 


5 


has  a  velocity  different  form  that  of  the  freestream  is 
allowed  to  pass  through  the  pressure  field  of  the  rotor 
blade  passage.  The  presence  of  the  freestream  fluid  is 
ignored.  Thus  no  mixing  or  interaction  between  the 
particle  and  the  freestream  is  considered.  However  with 
the  technique  developed  such  interactions  can  be  intro¬ 
duced  through  the  specification  of  the  velocity  and  the 
relative  flow  angle  evolution  of  the  particle.  Generally 
particles  that  have  a  velocity  defect  (as  apposed  to  a 
velocity  excess)  are  of  the  most  interest,  so  the  analy¬ 
sis  was  developed  in  terms  of  low  momentum  particles. 

The  analysis  does  however  apply  equally  well  to  particles 
that  have  a  velocity  greater  than  that  of  the  freestream. 
This  should  be  kept  in  mind  throughout  the  explanation  of 
the  calculation  technique.  The  particle  approach  helped 
a  great  deal  in  obtaining  an  intuitive  insight  into  what 
parameters  were  important  for  the  calculation  of  the  part¬ 
icle  trajectory.  The  author  relied  heavily  on  this  intui¬ 
tive  approach  to  lead  to  the  correct  mathamatical  formu¬ 
lation.  The  importance  of  correctly  modelling  the  physical 
behavior  of  the  particle  seemed  to  dictate  such  an  ap¬ 
proach,  making  mathamatical  derivations  and  justifications 
of  secondary  importance.  Thus  few  such  rigorous  matha¬ 
matical  principles  will  be  presented  in  this  analysis. 


2.  MECHANISMS  OF  ENERGY  TRANSFER 


The  fundamental  equation  for  the  energy  trans¬ 
fer  in  a  turbomachine  is  the  Euler  pump  equation 


AH=  VsVuZ~  lIYu,  f2-1' 


The  enthalpy  addition  in  a  compressor  is  dependent  upon 
the  change  in  tangential  velocity  from  inlet  to  exit  and 
on  the  uheel  speed.  Basically  two  different  mechanisms 
for  the  energy  transfer  exist.  They  are  the  residence 
time  of  the  particle  uithin  the  blade  rou  (resulting  in 
tangential  velocity  changes)  and  the  radial  shift  of  the 
particle  (resulting  in  changes  in  the  local  uheel  speea). 
Due  to  the  associated  axial  velocity  defect  of  a  lou 
momentum  particle  the  time  necessary  for  the  particle 
to  pass  through  the  rotor  blade  passage  is  greater  than 
that  of  the  freestream  fluid.  Smith  (REF  7)  shoued 
that  the  instantaneous  enthalpy  addition  for  an  inviscid 
fluid  particle  can  be  expressed  as 


DH  _  L 
Dt  ~  P  bt 


(2.2) 


where 

H  :  stagnation  enthalpy  per  unit  mass 
•£,  :  time 

p  :  density 
p  :  static  pressure 


7 


In  a  compressor  the  right  hand  side  of  this  equation  is 
positive,  so  the  energy  transfer  increases  as  the  resi¬ 
dence  time  for  a  particle  becomes  greater.  This  effect 
can  perhaps  be  seen  more  clearly  by  examining  the  velo¬ 
city  triangles  for  a  lou  momentum  particle  and  the  free- 
stream.  Figure  3a  is  such  a  velocity  triangle  for  a 
particle  having  a  defect  in  the  absolute  velocity  at  the 
enterance  of  a  rotor  blaae  row.  Note  that  the  lou  momen¬ 
tum  particle  has  a  smaller  absolute  tangential  velocity 
than  the  freestream  and  also  that  the  relative  flou  angle 
of  the  particle  is  larger  than  that  of  the  freestream. 

Thus  the  rotor  uill  "see"  such  lou  momentum  particles  as 
fluid  entering  the  blade  rou  uith  a  high  positve  incidence. 
It  is  assumed  throughout  this  analysis  that  a  particle 
regardless  of  the  magnitude  of  the  velocity  defect  or 
excess  uill  leave  the  rotor  blaae  rou  at  the  same  rela¬ 
tive  flou  angle  as  the  mean  freestream  flou.  The  rotor 
exit  velocity  triangle  for  the  freestream  and  a  lou 
momentum  particle  is  shoun  in  Figure  3b.  Here  it  can  be 
seen  that  due  to  associated  defect  in  the  axial  velocity 
the  lou  momentum  particle  has  a  larger  absolute  tangen¬ 
tial  velocity  than  the  freestream.  From  Eq.  2.1  it  can 
be  seen  that  (assuming  equal  uheel  speeds)  the  lou  momen¬ 
tum  particle  uill  have  a  larger  enthalpy  increase  than  the 
freestream.  As  the  magnitude  of  the  velocity  defect  of 
the  particle  becomes  greater  (i.e.  as  the  axial  velocity 
of  the  particle  is  reduced)  this  difference  in  the  energy 
transfer  for  the  particle  and  the  freestream  uill  become 


larger.  The  above  example  was  for  a  particle  entering  the 
blade  row  with  an  absolute  velocity  defect,  A  second  pos¬ 
sibility  exists  and  that  is  for  a  particle  entering  the 
blade  row  with  the  freestream  velocity  but  due  to  viscous 
(or  other)  effects  in  the  blade  passage  it  leaves  the 
blade  row  with  a  defect  in  the  exit  relative  velocity. 

The  inlet  and  exit  velocity  triangles  for  such  a  situa¬ 
tion  are  shown  in  Figure  4.  Again  due  to  the  velocity 
defect  the  particle  leaves  the  blade  row  with  a  larger 
absolute  tangential  velocity  and  thus  will  have  a  larger 
enthalpy  increase. 

Up  to  this  point  the  influence  of  changes  in 
the  wheel  speed  has  been  neglected.  From  Eq .  2.1  and 
the  velocity  triangles  of  Figure  3.  and  4.  the  influence 
of  the  wheel  speed  can  be  seen.  The  wheel  speed  increases 
proportionally  with  the  radial  location.  For  two  part¬ 
icles  with  equal  inlet  and  exit  tangential  velocities 
the  particle  with  the  larger  wheel  speed  will  have  the 
larger  enthalpy  increase.  Note  that  the  magnitude  of 
the  wheel  speed  also  influences  the  resulting  tangential 
velocity,  thus  having  a  double  effect  on  the  enthalpy 
increase  (i.e.  enthalpy  rise  increases  with  the  square  of 
the  wheel  speed).  The  influence  of  changes  in  wheel  speed 
on  the  energy  transfer  is  an  important  aspect  of  this  an¬ 
alysis  because  the  radial  pressure  field  of  the  blade  row 
will  tend  to  cause  low  momentum  fluid  particles  to  migrate 
radially  outward.  Throughout  this  analysis  it  is  vital 


to  keep  in  mind  the  relationship  of  the  velocity  compo¬ 
nents  and  the  uheel  speed  to  each  other  in  the  velocity 
triangles  for  the  particle  and  the  freestream. 


There  are  tuo  factors  that  will  determine  the 
magnitude  of  the  radial  migration  of  the  lou  momentum 
fluid  particles.  First,  the  magnitude  of  the  imbalance 
of  the  pressure  forces  and  the  centrifugal  forces  deter¬ 
mines  the  size  of  the  resultant  force  on  the  particle 
relative  to  the  freestream.  The  magnitude  of  this  force 
is  dependent  upon  the  difference  between  the  tangential 
velocity  of  the  particle  ana  the  tangential  velocity  of  the 
freestream  at  the  same  radial  location.  The  second  para¬ 
meter  influencing  the  amount  of  radial  shift  is  the  time 
this  net  force  is  applied  to  the  particle.  The  greater 
the  length  of  time  the  force  is  applied  the  larger  will 
be  the  deflection  of  the  particle  from  a  freestream 
streamline  trajectory.  This  is  the  so  called  resioence 
time  of  the  particle  and  it  is  directly  dependent  upon 
the  axial  velocity  of  the  particle.  Thus  the  size  of 
the  velocity  defect  effects  not  only  the  magnitude  of 
the  force  imbalance  of  the  particle  but  also  the  time 
this  force  is  applied.  The  relative  flow  angle  influences 
the  magnitude  of  the  tangential  velocity  and  thus  the  size 
of  the  applied  force  on  the  particle. 
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3.  THE  COMPUTER  PROGRAM 


Prior  to  investing  a  great  deal  of  time  and 
effort  creating  an  elaborate  computational  scheme  for 
calculating  the  trajectory  of  a  low  momentum  particle 
a  very  simple  computer  program  was  written  to  test  the 
feasibility  of  such  a  calculation.  It  was  also  possible 
to  gain  some  insight  into  the  assumptions  necessary  for 
and  the  restrictions  of  this  type  of  calculation  from 
an  initial  test  program. 


The  calculation  of  the  trajectory  of  a  particle 
is  centered  around  a  radial  equilibrium  equation.  For 
the  initial  program  the  simplest  form  of  the  radial 
equilibrium  equation  was  used 

±2l=&*  (3.-1) 

This  equation  is  a  statement  of  the  assumption  that  the 
centrifugal  forces  generated  by  the  freestream  fluid 
"circling"  the  compressor  axis  are  in  equilibrium  with 
the  radial  pressure  field.  Lou  momentum  particles  in 
the  flow  field  are  subjected  to  the  same  radial  pressure 
field,  but  due  to  their  lower  momentum  the  trajectory 
of  these  particles  will  be  influenced  more  by  the  radial 
pressure  field.  The  low  momentum  particles  will  "seek" 


a  radial  location  where  the  pressure  forces  and  centri¬ 
fugal  forces  are  again  balanced.  For  particles  with  an 


excess  of  tangential  velocity  (normally  associated  uith 
low  momentum  particles)  this  equilibrium  location  will 
occur  at  a  larger  radius.  This  equilibrium  equation  is 
the  only  velocity  related  equation  used  in  the  initial 
analysis  so  it  is  necessary  to  make  two  additional  assump¬ 
tions  in  order  to  solve  for  the  two  remaining  velocity 
components.  The  first  of  these  assumptions  is  that  the 
velocity  ratio,  (T  ,  which  is  the  ratio  of  the  particle 
axial  velocity  to  the  freestream  axial  velocity  is  a 
constant  through  the  entire  blade  row.  The  second 
assumption  necessary  concerns  the  evolution  of  the  rela¬ 
tive  flow  angle  of  the  particle.  In  this  initial  analysis 
only  defects  in  entering  absolute  velocity  were  considered 
From  Figure  3. a  it  can  be  seen  that  the  inlet  relative 
flow  angle,  fi  ,  is  established  by  the  freestream  flow 
direction  and  the  velocity  ratio,  (T  .  The  resulting 
relative  flow  angle  for  the  particle  can  be  consider¬ 
ably  larger  than  the  freestream  angle  at  the  inlet,  but 
these  two  angles  are  assumed  to  be  equal  at  the  blade 
exit.  Thus  it  is  necessary  to  specify  how  the  relative 
flow  angle  of  the  particle  evolves  through  the  blade  row. 

A  simple  linear  variation  was  used.  That  is,  when  the 
particle  is  one-half  of  the  way  through  the  blade  row 
the  difference  between  the  relative  flow  angles  of  the 
particle  and  the  freestream  is  assumed  to  be  one-half  of 
what  it  was  at  the  inlet  of  the  blade  row.  These  two 
assumptions  were  made  to  keep  the  initial  test  program 
as  simple  as  possible.  Uith  the  equilibrium  equation  of 
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3.1  and  these  assumptions  the  initial  program  uas  written 
and  tested.  The  freestream  velocity  field  uas  generated 
by  a  Katsanis-l'lc.Mally  through-flow  calculation  and  used 
as  input  data  for  the  program.  The  results  of  this  simpl¬ 
ified  approach  were  encouraging  and  enlightening.  There 
were  some  discovered  limitations  of  the  simple  approach 
but  the  technique  seemed  to  warrant  further  improvements 
and  refinements.  These  improvements  would  involve  the 
equilibrium  equation  and  the  limitations  of  the  velocity 
and  relatve  flow  angle  evolution  used. 

The  final  program  consists  of  a  MAIN  deck  and 
four  supplimentary  subroutines.  All  of  the  trajectory 
and  energy  calculation  are  done  in  the  MAIN  deck.  Two 
interpolation  subroutines  are  used  to  obtain  velocity/ 
space  interpolations  (LININT)  and  interpolations  of  var¬ 
ious  parameters  that  are  a  function  of  one  variable  only 
(INTERPOL).  The  two  remaining  subroutines  are  concerned 
with  a  boundary  layer  thickness  calculation  (BONLAY)  and 
the  resulting  boundary  layer  velocity  profile  (UELRmTIO). 
Because  the  principle  objective  of  this  project  was  the 
development  of  the  calculation  technique, the  procedure 
will  be  explained  in  some  detail.  For  persons  interested 


only  in  using  the  program  the  input  and  output  are  ex¬ 
plained  in  Appendix  I  and  II  along  with  a  program  list¬ 
ing  and  examples  of  input  and  output  data. 


1  3 


3.1  PIAIN  Deck 

The  results  of  the  initial  program  indicated 
a  serious  limitation  of  the  equilibrium  equation  used. 

Uhen  a  particle  uas  introduced  uith  a  velocity  ratio  equal 
to  unity  (i.e.  Q~  =  1.00)  the  trajectory  calculated  uas 
that  of  a  straight  line  in  the  meridional  plane.  The 
slope  of  this  trajectory  uas  determined  by  the  initial 
radial  velocity  uhich  then  remained  constant  through  the 
blade  rou.  The  equilibrium  equation  did  not  take  into 
account  the  radial  acceleration  of  the  freestream  (i.e. 
changes  in  radial  velocity  along  a  freestream  stream¬ 
line).  This  simple  equilibrium  equation  took  into  account 
only  the  acceleration  of  the  particle  relative  to  the 
freestream  (uhich  is  zero  for  ( T  =  1.00).  The  reason  for 
this  result  can  be  seen  by  examining  the  radial  component 
of  the  Euler  Equations 

Fhet  _  1  ,  _Vu  _  DVr  n\ 

e  ~  9  R  D± 

This  equation  is  Neuton's  2nd  Lau  applied  to  a  fluid 
particle,  equating  the  sum  of  the  forces  (per  unit  mass) 
to  the  net  acceleration  of  the  particle.  It  can  be  seen 
that  use  of  the  initial  equation  for  the  radial  equilibrium 
(3.l)  assumed  that  no  net  acceleration  of  the  freestream 
existed.  This  is  not  the  case  if  the  radial  velocity  of 
the  freestream  varies  along  a  streamline  (due,  for  example, 
to  meridional  curvatures). 
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It  uas  also  realized  that  for  a  velocity  ratio 
of  unity  (  (T  =  1.00)  the  calculation  scheme  must  be  able 
to  properly  calculate  the  radial  location  of  a  freestream 
streamline  through  the  blade  rou.  Clearly  3.1  is  not 
capable  of  doing  this  because  there  were  large  radial  var¬ 
iations  in  the  freestream  due  mainly  to  the  meridional 
curvatures.  Thus  a  first  "test  case"  for  any  neu  calcula¬ 
tion  scheme  will  be  "can  it  properly  calculate  the  rauial 
location  of  a  freestream  streamline  through  the  blade 
rou"  . 


H  brief  explanation  of  the  equilibrium  equation 
used  in  the  final  program  will  now  be  given.  The  radial 
equilibrium  equation  of  3.2  must  apply  for  all  particles 
in  the  flou  field.  Writing  the  equation  for  a  freestream 
(FS)  particle  and  a  lou  momentum  (p)  particle  yields 

DVk  \  _  VC* _ i  j>P  ^  . 

ot  R  Is  P  aR  ^  3,3 


QVfi  \  -  Vo.  ) _ i_  _bP  \ 

Oh  R  Jp  ?  bR  'P 


(3.4) 


The  radial  pressure  field  of  the  freestream  is  imposed 
upon  the  lou  momentum  particles.  Therefore, 


j_  aP  \ _ LIE) 

e  zr  >  p  « 


(3.5) 
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Rewriting  the  freestream  equilibrium  equation  as 

--L&)  =££)  -A)  (3.6) 

?  'fs  D-l'fs  R  'fs 

and  substituting  in  for  the  pressure  gradient  term  in  the 
equilibrium  equation  of  the  low  momentum  particle  yields 

ok.)  =  -b  y±.)  -  (3.7) 

O-t  Jp  Oi.  JFS  R  'p  R'FS 

From  this  equation  it  can  be  seen  that  the  total  radial 
acceleration  of  the  particle  is  equal  to  that  of  the  free¬ 
stream  plus  the  acceleration  of  the  particle  relative  to 
the  freestream  due  to  an  imbalance  of  the  centrifugal  ana 
pressure  forces.  If  the  centrifugal  terms  are  of  equal 
magnitude  for  the  particle  and  freestream,  the  net  raaial 
acceleration  of  the  particle  is  equal  to  that  of  the  free¬ 
stream  and  thus  should  follow  the  same  trajectory  (or 
streamline)  as  a  freestream  particle. 

It  was  found  that  a  freestream  meridional 
streamline  could  be  satisfactorily  approximated  by  a 
Taylor  series  expansion  for  the  radius,  R.  That  is, 

2  2 

R(x^o=  +  +  sL B)  AL  (3.8) 

cJ~£-  *  (j  -£.2  'j:  2. 

where  "I"  and  "1+1"  indicate  axial  computational  stations. 
This  equation  can  be  simplified  and  solved  by  recognizing 
that  the  first  derivative  term  of  Equation  3.B  is  the 
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radial  velocity  and  the  second  derivative  term  is  the 
radial  accaleration .  Approximating  the  acceleration 

term  by  a  forward  difference  scheme  and  neglecting  all 

2 

terms  of  the  )  and  smaller  yields, 

R(X-H)  =  +  ^(X)A-k  +[[£ AL  (3.9) 

*  L  At  j  2 

which  can  be  further  simplified  to 

R(**0S  R/i)  +  VS(-r*0+VR(x)J  (3>10) 

This  trajectory  equation  states  that  the  new  radial  lo¬ 
cation  (i.e.  at  the  new  axial  station)  is  equal  to  the 
old  radial  location  plus  the  average  radial  velocity  of 
the  particle  over  the  interval  times  the  time  the  part¬ 
icle  is  within  this  interval.  The  resulting  trajectory 
for  a  freestream  particle  (i.e.  (J"  =  1.U0)  is  in  very 
good  agreement  with  the  radial  location  of  the  meridional 
streamline  calculated  by  the  through-flow  calculation. 
Solving  equation  3.10  involves  an  iterative  solution  be¬ 
cause  Vpj(l  +  l)  is  a  function  of  the  new  radial  location, 

R(  I  +  'l  )  . 


The  equilibrium  equation  is  again  the  only 
equation  that  relates  the  trajectory  of  the  particle 
to  the  freestream  flow  field.  Thus  it  is  possible  to 
calculate  only  one  new  velocity  component.  Auxilary 
equations  or  assumptions  are  necessary  for  the  two 
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remaining  velocity  components.  As  with  the  initial  pro¬ 
gram  it  uas  decided  to  specify  as  input  the  additional 
pieces  of  information  but  remove  some  of  the  limiting 
restrictions  of  the  initial  program.  The  two  additional 
parameters  are  again  the  evolution  of  the  velocity  ratio, 
c r  ,  and  the  relative  flou  angle  of  the  particle. 

Every  effort  uas  made  to  keep  this  program  as  flexible 
as  possible  by  giving  the  user  complete  control  over 
these  parameters  via  an  easily  changed  input.  The  values 
of  the  velocity  ratio  and  relative  flou  angle  can  be 
specified  as  arbitrary  functions  of  the  axial  distance 
through  the  blade  rou,  or  for  relative  velocity  defects 
(such  as  surface  velocity  distributions  or  boundary 
layer  particles)  calculated  within  the  program  from  the 
surface  velocity  data  or  the  boundary  layer  model. 


The  equilibrium  equation  of  3.7  is  solved 
using  a  forward  difference  scheme  yielding 


2  2 

L  At  -lp  L  JZ  JFS  R  ?  r'fs 


3.11) 


A  physical  interpretation  of  this  equation  is  easier 
when  equation  3.11  is  written  as 


^  -  *«>  ♦£* (IH>  - « ix)ls  li  i  T V  hh 


(3.12) 


The  new  radial  velocity  (at  the  next  axial  station) 

V  (i+l)  is  equal  to  the  old  radial  velocity  plus  the 

n 


18 


change  in  freestream  radial  velocity  "weighted"  by  the 
ratio  of  residence  time  for  a  low  momentum  particle  and 
the  freestream  plus  an  aoditional  term  due  to  the  dif¬ 
ference  in  centrifugal  and  pressure  forces.  From  this 
equation  it  is  easy  to  see  the  influence  of  the  residence 
time  and  the  tangential  velocity  difference  on  the  rate 
of  change  of  the  radial  velocity  of  the  particle. 

Equation  8.12  ana  3.10  are  used  to  calculate 
the  velocity  and  radial  location  of  the  particle  at  a 
new  axial  calculation  station.  From  the  known  freestream 
velocity  values  the  axial  velocity,  \l ^  ,  and  the  tangen¬ 
tial  velocity,  V ,  of  the  particle  can  be  determined 
using  the  supplied  input  data  for  the  velocity  ratio  ana 
relative  flow  angle  evolution.  However,  because  these 
values  are  functions  of  the  new  radial  location  it  is 
necessary  again  to  "guess"  a  new  radial  location  and 
solve  equation  3.12  ans  3.10  iteratively  until  the 
guessed  radial  location  and  the  calculated  value  are 
within  a  specified  tolerance. 

Once  the  complete  trajectory  of  the  particle 
is  obtained  the  energy  transfer  calculations  can  be  made. 
These  calculations  are  made  for  the  three  streamlines 
shown  in  Figure  5.  Streamline  i\  is  the  freestream  stream¬ 
line  coinciding  with  the  radial  location  of  the  particle 
at  the  inlet  of  the  blade  row.  Streamline  B  is  the  free¬ 
stream  streamline  coinciding  with  the  radial  location  of 


the  particle  at  the  exit  of  the  blade.  The  temperature 
difference  between  streamline  A  and  the  particle  repre¬ 
sents  the  additional  energy  transfer  due  to  the  different 
velocity  evolution  of  the  particle.  The  temperature  dif¬ 
ference  between  streamline  B  and  the  particle  indicates 
the  size  of  the  temperature  variations  in  the  flow  field 
that  will  exist  at  the  exit  of  the  rotor. 

3 . 2  Subroutine  LINIIMT 

This  subroutine  is  a  three-dimensional  inter¬ 
polation  routine  used  to  obtain  velocity  values  of  the 
freestream  at  a  specified  meridional  location,  from  the 
streamline  values  obtained  from  the  through- flow  calcu¬ 
lation.  The  through-flow  calculation  provides  freestream 
velocity  values  along  seven  meridional  streamlines  (hub, 
10$,  30$,  50$,  70$,  90$,  and  tip)  with  10  to  15  values 
within  the  blade  row  for  each  streamline.  This  routine 
was  obtain  from  an  existing  program. 

3 . 3  Subroutine  INTERPOL 

This  subroutine  is  also  an  interpolation  rou¬ 
tine.  It  is  used  to  interpolate  between  values  of  surface 
velocity,  relative  flow  angle,  and  velocity  ratio  that 
are  supplied  as  functions  of  the  axial  distance  only. 

This  interpolation  routine  was  written  by  the  author. 
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3 . 4  Subroutine  BO  NLA  Y 


This  routine  is  a  simple  boundary  layer  cal¬ 
culation  used  to  obtain  the  total  boundary  layer  thick¬ 
ness  evolution,  the  skin  friction  coefficient,  and  a 
velocity  profile  parameter  (Lefoll's  L)  for  the  rotor 
blade.  Emperical  correlations  are  used  for  the  conven¬ 
tional  shape  factor  (H,^)  anc>  sl<^n  friction  coef¬ 

ficient  (REF  11 ).  An  additional  emperical  correlation 
uas  used  (see  REF  12)  to  obtain  the  total  boundary  layer 
thickness  from  the  calculated  momentum  thickness.  This 
value  is  dependent  upon  the  assumed  boundary  layer  velo¬ 
city  profile  (Lefoll).  Some  difficulty  arose  in  defin¬ 
ing  a  starting  condition  for  the  calculation.  A  simple 
fundamental  equation  uas  derived  and  used  to  establish 
a  flat  plate  type  starting  condition.  The  boundary 
layer  thickness  calculation  uas  a  finite  differnece 
scheme  (using  a  central  difference)  for  the  intergral 
boundary  layer  equation.  The  calculation  procedure  as¬ 
sumed  the  boundary  layer  to  be  turbulent  and  no  condition 
uas  used  to  detect  separation.  Thus  care  must  be  used 
uhen  calculating  boundary  layers  uhere  separation  may 
occur . 


3 . 5  Subroutine  VELRATIQ 


Once  the  total  boundary  layer  thickness  evolu¬ 


tion  is  knoun  it  is  possible  to  calculate  the  velocity 
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defect  of  a  boundary  layer  particle  at  any  location  with- 
in  the  boundary  layer  for  an  assumed  velocity  profile. 

To  keep  the  boundary  analysis  as  simple  as  possible 
it  was  decided  to  use  the  Lefoll  (REF  lli)  boundary  layer 
profile.  This  velocity  profile  can  be  represented  by 

2 

+  L(27)~  371*+  I )J  (3.13) 

From  a  specified  location  within  the  boundary  layer  it 
is  possible  to  obtain  the  local  velocity  ratio  for  a 
particle  at  that  location  in  the  boundary  layer.  An 
assumptions  was  necessary  concerning  the  shape  of  the 
streamline  of  a  boundary  layer  particle  in  the  blade- 
to-fclade  plane.  The  assumption  made  was  that  a  boundary 
layer  particle  will  travel  on  a  streamline  that  is  a 
constant  distance,  y,  from  the  surface  of  the  blade. 

This  distance  is  spcified  as  a  fraction  of  the  exit 
boundary  layer  thickness,  ^  .  Thus  it  is  possible  to 
choose  how  "deep"  in  the  boundary  layer  the  particle  of 
interest  is.  As  shown  in  Figure  6.,  for  y/g£  =  1.00  the 
velocity  of  the  particle  will  be  equal  to  the  freestream 
surface  velocity.  As  y/<££  becomes  smaller  the  particle 
enters  the  boundary  layer  progressively  nearer  the  blade 
leading  edge  and  thus  the  particle  will  leave  the  blade 
row  with  a  larger  velocity  defect.  The  velocity  defect 
is  not  only  a  function  of  the  specified  y/£g  but  also 
of  the  evolution  of  the  boundary  layer  thickness. 
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3 . 6  Modes  of  Operation 

The  program  uas  developed  to  handle  too  dis¬ 
tinctly  different  applications.  The  specification  of 
the  velocity  ratio  has  been  separated  for  absolute  velo¬ 
city  ratios  and  relative  velocity  ratios  (see  Figure  3. 
and  Figure  4.).  For  a  defect  in  the  absolute  velocity 
the  particle  enters  the  blade  rou  with  a  relative  flou 
angle  different  from  that  of  the  freestream.  This  rela¬ 
tive  flou  angle  difference  does  influence  the  trajectory 
of  the  particle  because  it  changes  the  tangential  velocity 
term  in  equation  3.12.  For  this  mode  of  operation  it  is 
necessary  to  supply  as  input  hou  the  relative  flou  angle 
progresses  from  the  inlet  value  to  the  exit  value.  In 
this  mode  of  operation  it  is  possible  to  model  inlet  flou 
nonunifarmities  such  as  casing  boundary  layers  and  up¬ 
stream  blade  uakes  entering  the  blade  rou  and  determine 
the  influence  of  these  flous  on  the  energy  transfer. 

Due  to  a  lack  of  time  and  the  additional  complexity  of 
establishing  uake  or  boundary  layer  decay  models  this 
mode  of  operation  uas  not  used  during  the  course  of  this 
project.  It  should  be  noted  houever  that  such  entrain¬ 
ment,  uake  decay  and  mixing  models  do  exist  in  the  liter¬ 
ature  and  could  be  incorporated  in  this  calculation 
procedure.  The  input  format  nou  is  such  as  to  allou  the 
user  to  quickly  and  easily  change  the  velocity  ratio  and 
flou  angle  evolution  data. 

1 

1 


The  second  mode  of  operation  is  for  velocity 


ratios  in  the  relative  velocity.  Here  the  relative  flow 
direction  is  not  influenced  by  the  velocity  ratio,  (T  . 

The  particle  will  have  the  freestream  relative  flow  angle 
everywhere  in  the  blade  row  making  the  specification  of 
the  relative  flow  angle  evolution  unnecessary.  One  fur¬ 
ther  simplification  was  made  for  this  mode  of  operation. 
All  blade- toblade  streamlines  (suction  surface  to  pressure 
surface)  are  considered  to  be  the  same  shape  as  the  mean 
streamline.  That  is,  the  relative  flow  angle  evolution 
for  all  blade- to-blade  streamlines  is  assumed  equivalent 
to  that  of  the  mean  freestream  which  is  essentially  the 
blade  camber  line.  Thus  the  blades  are  replaced  by  their 
camber  lines  only  as  shown  in  Figure  7.  For  thin  blades 
away  from  the  leading  edge  this  is  a  reasonable  approxi¬ 
mation.  The  through-flow  calculation  scheme  included  a 
blade- to-blade  calculation  as  well  to  obtain  the  blade 
surface  velocity  distributions.  Because  the  relative 
flow  angle  is  assumed  equal  for  all  blade- to-blade  stream¬ 
lines  the  trajectory  of  a  particle  becomes  dependent  only 
on  the  relative  velocity  ratio,  0"*  .  This  mode  of  opera¬ 
tion  is  useful  for  determining  the  influence  of  the  blaoe- 
to-blade  velocity  gradient  and  the  rotor  blade  boundary 
layers  on  the  energy  transfer  of  the  blade  row.  The 
freestream  (out  of  the  boundary  layers)  blade-to-blade 
streamline  of  interest  can  be  specified  in  the  input  data 
(Y/PITCH  can  range  from  0.0  on  the  suction  surface  to 
1.0  on  the  pressure  surface  with  the  mean  streamline 
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being  0.5).  An  interpolation  routine  is  used  to  inter¬ 
polate  between  the  mean  streamline  velocities  and  the 
surface  welocities.  This  mode  of  operation  was  used  ex¬ 
tensively  during  the  project  to  obtain  freestream  and 
boundary  layer  particle  trajectories  and  the  resulting 
rotor  exit  total  temperature  field  variations  for  a 


transonic  rotor 
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4.  RESULTS 

The  flou  field  of  a  transonic  rotor  uas  used 
to  demostrate  the  use  and  capability  of  the  program 
developed.  The  program  uas  used  to  calculate  the  in¬ 
fluence  of  the  b lad e- to-blad e  velocity  gradient  and  the 
boundary  layer  migration  on  the  rotor  exit  flou  field. 
The  results  uill  first  be  presented  for  the  freestream 
(i.e.  fluid  outside  the  blaae  boundary  layers)  surface 
particles  and  thenfor  the  blade  boundary  layer  fluid. 
Results  uill  be  shoun  at  span  locations  of  IQ70,  5Q>, 
and  90%  (based  of  flou  area).  The  calculated  trajec¬ 
tories  for  these  particles  and  the  influence  of  these 
trajectories  on  the  energy  transfer  of  the  rotor  uill 
then  be  shoun.  All  radial  shift  and  energy  transfer 
comparisons  are  made  uith  respect  to  the  mean  blade- 
to-blade  streamline  of  the  through-flou  calculation. 

The  analysis  uas  made  for  clean  inlet  flou  conditions 
only.  l\lo  attempt  uas  made  to  introduce  inlet  flou 
field  nonuniformities.  Finally,  a  short  summary  of 
hou  these  results  compare  uith  existing  test  data  uill 
be  given. 


4.1  Freestream  Surface  Velocity  Evolution 

Due  to  the  curvature  of  the  blade  channel  the 
surface  velocities  of  the  blade  uill  be  different  than 
the  mean  blade- to-blade  streamline  velocity.  The  mag- 
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nitude  of  these  velocity  differences  depends  upon  the 
blade  design  and  loading.  In  this  analysis  the  axial 
velocity  and  tangential  velocity  are  obtain  from  the 
freestream  velocity  values  (generated  by  the  through- 
flou  calculation)  and  the  boundary  layer  calculation. 

Only  the  radial  velocity  component  is  actually  calcu¬ 
lated  within  the  program.  In  Figure  8. a  the  axial 
velocity  evolution  for  the  mean,  suction  surface,  and 
pressure  surface  streamlines  is  shown  for  the  50%  span 
location.  The  velocity  distributions  are  shown  versus 
a  nondimensional  axial  distance,  0/a  representing  the 
blade  leading  edge  and  1 00;6  the  blade  trailing  edge. 

The  suction  surface  velocity  is  everywhere  greater  than 
that  of  the  mean  streamline  velocity  producing  a  velo¬ 
city  ratio,  ( T  ,  larger  than  one.  The  opposite  is  true 
of  the  pressure  surface  velocity,  is  everywhere  less 

than  one.  The  residence  time  variations  of  mean,  suc¬ 
tion  surface,  and  pressure  surface  particles  are  a  direct 
result  of  these  axial  velocity  differences.  From  the 
velocity  triangles  it  can  be  seen  that  (recall  for  this 
analysis  the  relative  flow  angles  of  all  blade- to-blade 
streamlines  are  assumed  equal  to  the  mean  streamline 
flow  angle)  the  suction  surface  particles  will  have  a 
smaller  absolute  tangential  velocity  and  pressure  sur¬ 
face  particles  a  larger  absolute  tangential  velocity 
than  the  mean  streamline.  The  absolute  tangential  velo¬ 
city  evolution  through  the  blade  row  for  each  of  these 
streamlines  is  shown  in  Figure  8.b,  again  for  50/6  span. 
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Because  the  differences  in  tangential  velocity  is  an 
indication  of  the  magnitude  of  the  force  imbalance  on  the 
particle  it  is  possible  to  obtain  an  idea  of  the  evolu¬ 
tion  of  this  force  on  the  particle  from  this  figure. 

The  larger  tangential  velocity  for  pressure  surface 
particles  will  tend  (see  eq .  3.12)  to  produce  a  larger 
radial  acceleration  of  the  particle,  uhile  the  opposite 
trend  occurs  for  suction  surface  particles.  The  result¬ 
ing  calculated  radial  velocity  evolution  is  shoun  in 
Figure  8.c.  As  expected  due  to  the  larger  residence 
time  and  tangential  velocity  the  pressure  surface  radial 
velocity  is  larger  than  that  of  the  mean  streamline.  The 
suction  surface  velocity  is  smaller  than  the  mean.  A 
few  remarks  can  be  made  concerning  these  radial  velocity 
evolutions.  From  1  0%  to  40>t  meridional  distance  the  tan¬ 
gential  velocity  difference  betueen  the  pressure  surface 
and  the  mean  streamline  is  the  greatest.  The  resulting 
rate  of  change  of  the  pressure  surface  radial  velocity 
is  also  the  largest.  From  40%  to  70>S  meridional  dis¬ 
tance  the  tangential  velocity  difference  (indication  of 
force  imbalance)  becomes  smaller,  but  still  non-zero. 
However  the  radial  velocity  of  the  pressure  surface 
particle  actually  decreases  in  this  range.  This  can  be 
explained  by  noting  the  rather  large  rate  of  change  of 
the  radial  velocity  for  the  mean  streamline.  Hecall, 
the  change  in  radial  velocity  of  a  particle  is  related 
to  the  change  of  radial  velocity  of  the  freestream 
"weighted"  by  the  ratio  of  resiaence  times.  Here  the 
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large  rate  of  change  of  the  freestream  radial  velocity 
and  the  increased  residence  time  of  the  pressure  surface 
particle  combine  to  overcome  the  centrifugal  force  im¬ 
balance  that  trys  to  increase  the  particles  radial  velo¬ 
city.  This  is  essentially  the  influence  of  streamline 
curvature  on  the  low  momentum  particle.  For  the  suction 
surface  particles  all  tendencies  are  to  reduce  its  radial 
velocity  relative  to  the  mean  freestream. 

Similar  axial,  tangential,  and  resulting  cal¬ 
culated  radial  velocity  evolutions  are  shown  in  Figure  9. 
and  Figure  10.  for  the  10 %  (near  the  hub)  span  location 
and  the  90%  (near  the  tip)  span  location.  Figure  9.c 
shows  clearly  the  influence  of  the  high  streamline  curva¬ 
ture  near  the  hub  on  lower  and  higher  momentum  particles. 
Again  the  centrifugal  force  imbalance  tending  to  push 
the  pressure  surface  particles  upward  (thus  increasing 
the  radial  velocity)  is  overcome  by  the  strong  change  in 
the  mean  streamline  radial  velocity.  The  suction  surface 
particles  because  of  their  greater  through  flow  velocity 
are  less  effected  by  the  high  streamline  curvature  and 
thus  have  a  larger  radial  velocity  than  the  mean  stream¬ 
line  through  the  entire  blade  row.  The  higher  wheel 
speed  at  the  tip  (Figure  10.)  tends  to  produce  larger 
tangential  velocity  differences  and  thus  larger  differ¬ 
ences  between  the  radial  velocity  for  pressure  surface, 
mean,  and  suction  surface  particles.  The  radial  velo¬ 
city  of  pressure  surface  particles  is  beginning  to  become 
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significant  at  the  tip. 

4  . 2  Boundary  Layer  Velocity  Evolution 

The  boundary  layer  results  are  presented  for 
four  different  boundary  layer  locations.  These  locations 
are  specified  as  a  ratio  of  the  distance  from  the  blade 
surface,  y,  to  the  exit  boundary  layer  thickness,  $£  . 

A  value  of  1 .00  corresponds  to  the  edge  of  the  boundary 
layer  region  and  thus  is  identical  to  the  surface  velo¬ 
city  evolutions  presented  above.  Values  of  y/£^  of 
0.25,  0.05,  and  0.01  indicate  streamlines  successively 
closer  to  the  blade  surface  ana  thus  "deeper"  within  the 
boundary  layer. 

The  influence  of  the  boundary  layer  lo¬ 
cation  on  the  axial  velocity  evolution  for  the  suction 
surface  is  shown  in  Figure  11. a  for  50>j  span.  Note  that 
over  the  first  50>£  axial  distance  of  the  blade  most  of 
the  boundary  layer  fluio  still  has  an  excess  of  axial 
velocity  relative  to  the  mean  streamline.  The  tangen¬ 
tial  velocity  evolution  for  these  particles  is  shown  in 
Figure  11. b,  and  the  calculated  radial  velocity  evolution 
in  Figure  11. c.  As  expected  the  radial  velocities  of 
the  suction  surface  boundary  layer  particles  are  smaller 
than  the  mean  streamline  radial  velocity  over  most  of  the 


blade 
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similar  pressure  surface  boundary  layer  data 
are  shoun  in  Figure  12.  Again  results  are  shoun  for  four 
different  boundary  layer  locations.  The  freestream  sur¬ 
face  velocity  of  the  pressure  surface  is  already  con¬ 
siderably  less  than  that  of  the  mean  stramline.  As 
particles  deeper  uithin  the  boundary  layer  are  considered 
this  velocity  defect  becomes  even  larger.  From  the 
velocity  triangles  for  the  boundary  layer  particles  the 
tangential  velocity  evolution  is  obtained  and  is  shoun 
in  Figure  12. b.  Recall,  the  axial  velocity  ratio  is  an 
indication  of  the  residence  time  and  the  tangential  velo¬ 
city  excess  is  an  indication  of  the  radially  outuard 
force  on  the  particle.  The  particle  deepest  uithin  the 
boundary  layer  have  not  only  the  largest  residence  time 
but  also  the  greatest  force  imbalance.  Thus  the  high 
radial  velocities  shoun  in  Figure  12. c  for  these  particles 
is  not  too  surprizing.  The  largest  calculated  radial 
velocity  in  the  boundary  layer  is  approximately  1 0$  of 
the  freestream  axial  velocity. 

4 . 3  Particle  Trajectory  and  Radial  Shift 

The  meridional  trajectory  for  the  freestream 
surface  particles  and  the  boundary  layer  fluid  uill  be 
presented  in  this  section.  The  radial  location  of  the 
particle  is  dependent  not  only  on  the  evolution  of  the 
radial  velocity  of  the  particle  through  the  blade  but 
also  on  the  time  the  particle  remains  uithin  the  blade 
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row.  The  trajectories  for  freestream  suction  surface 
and  pressure  surface  particles  are  shown  in  Figure  13. 
for  the  1 0;i ,  50%,  and  90;‘!  span  locations.  The  twisting 
of  the  strearashest  as  it  progress  through  the  blade  rou 
can  be  seen.  The  amount  of  twist  varies  from  about  1  mm 
at  the  hub  to  2.5  mm  at  the  tip  with  the  nressure  surface 
always  higher  than  the  suction  surface. 

The  meridional  trajectories  for  suction  surface 
boundary  layer  particles  are  shown  in  Figure  14.,  also 
for  the  10%,  50% ,  and  90%*  span  locations.  There  is  very 
little  variation  in  the  particle  paths  for  the  suction 
surface  boundary  layer  fluio.  The  racial  shift  is  small 
and  nearly  constant  over  the  span. 

Similar  trajectory  results  are  shown  in  Figure 
15.  for  the  pressure  surface  boundary  layer  particles. 
Because  of  the  larger  force  imbalance  ano  the  greater 
residence  time  there  is  a  significant  radial  migration 
of  the  pressure  surface  boundary  layer  fluid.  The  mag¬ 
nitude  of  this  radial  shift  varies  over  the  blade  span, 
from  about  4  mm  at  the  hub  to  over  12  mm  at  the  tip. 

The  deepest  boundary  layer  particle  considered  ( y/Sg  = 
0.01 )  reached  the  shroud  casing  radius  before  leaving 


the  blade  rou. 


4.4  L  ne  rn  y  T  rans  f e  r 


ns  mentioned  in  Chapter  2  the  radial  shift 
can  play  an  important  role  in  determining  the  magnitude 
of  the  energy  transfer  to  a  particle.  In  Figure  16.  a 
circumferential  i/ieu  of  the  blade- to-blade  flow  field 
is  shown.  The  circumferential  distance  is  broken  into 
three  small  segments  showing  the  suction  surface  bound¬ 
ary  layer,  the  mean  streamline,  and  the  pressure  surface 
boundary  layer.  Zero,  on  the  horizontal  axis,  corresponds 
to  the  suction  surface  of  the  blade  ano  10 to  the 
pressure  surface.  Note  the  different  scale  used  for  the 
suction  surface  and  pressure  surface  boundary  layers. 

From  this  view  the  twisting  of  the  blaae-to-blaae  stream - 
sheet  can  clearly  be  seen.  Also  the  radial  migration  of 
pressure  surface  and  suction  surface  boundary  layer  part¬ 
icles  is  shown  in  the  figure.  The  bar  graphs  represent 
the  temperature  increase  these  fluid  particles  recieve 
in  addition  to  that  of  the  mean  streamline  due  to  their 
different  velocity  evolution.  The  overall  mean  free- 
stream  temperature  rise  through  the  blade  row  is  approx¬ 
imately  22°  C. 


Figure  16. a  shows  the  circumferential  varia¬ 
tions  in  the  flow  field  at  50, ^  span.  It  can  be  seen  that 
freestream  suction  surface  fluid  recieves  the  smallest 
energy  addition  while  pressure  surface  boundary  layer 
fluid  recieves  the  largest.  The  temperature  variation 
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is  shown  to  be  4°  -  5°  C.  at  this  span  location.  figure 
16. b  and  figure  16. c  show  similar  data  for  1  0;4  span  ana 
90;b  span.  Due  to  the  smaller  radial  shifts  ana  lower 
wheel  speed  at  1Q>  span  (figure  16. b)  the  circumferential 
temperature  variations  are  very  small.  At  90/j  span  how- 
ever  the  large  radial  shifts  combine  with  the  high  wheel 
speed  to  produce  the  largest  circumferential  temperature 
variations  in  the  blade  row. 

4.5  Relation  to  Test  Data 


Unfortunately  no  test  data  are  available  for 
the  transonic  rotor  used  in  the  analysis  above.  There 
are  however  some  experimental  results  available  for 
similar  machines.  High  speed  data  are  necessary  for 
temperature  variation  comparisons,  but  rotor  wake  raaiax 
velocity  comparisons  can  be  made  for  both  low  speed  ana 
high  speed  rotor  data. 

Several  experimenters  (REf  2,3)  have  reported 
measuring  rotor  wake  radial  velocities  on  the  order  of 
10;»  to  1  5;a  of  the  axial  velocity  downstream  of  a  rotor. 
This  is  the  order  of  magnitude  of  radial  velocity  ob¬ 
tained  for  particles  deep  within  the  pressure  surface 
boundary  layer. 

Leboeuf  (REf  9)  ana  Kool  (REf  10)  made  circum¬ 
ferential  temperature  field  measurements  behind  a  high 


speed  rotor  at  VKI.  They  reported  temperature  variations 
from  the  freestream  flow  to  the  blade  wake  fluid  to  be 
2°  -  3°  C  at  the  blade  mid-span.  The  rotor  used  for 
these  measurements  uas  not  as  highly  loaded  as  the  rotor 
used  in  the  above  analysis.  Kerrebrock  and  ('likola  jczak 
(REF  b)  report  circumferential  temperature  variations 
at  the  exit  of  a  stator  rou  dounstream  of  a  transonic 
rotor  to  be  on  the  order  of  "1  Q%  of  the  overall  stage 
temperature  rise.  Again  the  calculated  values  obtained 
from  the  program  compare  uell,  on  an  order  of  magnitude 
basis,  with  these  experimental  values. 

Although  it  is  difficult  to  drau  to  many  con¬ 
clusions  from  this  limited  comparison  it  is  felt  that 
the  calculation  technique  developed  above  does  produce 
results  that  are  realistic  and  are  similar  to  experi¬ 
mentally  determined  values. 
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5.  CONCLUS  IlJNS 


The  derivation  and  development  of  a  technique 
to  calculate  the  trajectory  of  lou  momentum  fluid  particles 
and  their  influence  on  the  energy  transfer  o^  a  rotating 
blade  rou  has  been  presented.  The  resulting  values  of 
radial  velocity  for  the  boundary  layer  fluid  and  the  as¬ 
sociated  increased  temperature  rise  for  these  particles 
agree  well  with  available  test  data. 

A  tuisting  or  uarpping  of  a  blade- to-blade 
streamsheet  through  a  rotor  blade  rou  uas  shoun  to  occur 
due  to  the  velocity  gradients  in  the  blade  channel.  The 
suction  surface  freestream  fluid  tends  to  go  aounuard 
uhile  the  pressure  surface  fluid  goes  radially  outuaro. 

A  quantitative  value  for  this  streamsheet  tuisting  uas 
given  for  a  transonic  rotor.  The  amount  of  tuist  is 
greatest  at  the  blade  tip. 

The  radial  migration  of  rotor  blade  boundary 
layer  fluid  uas  shoun  to  occur.  The  radirl  shift  of  this 
fluid  is  more  significant  in  the  pressure  surface  bouno- 
ary  layer  than  in  the  suction  surface  boundary  layer. 

The  migration  of  the  boundary  layer  fluid  uas  shoun  to 
be  influenced  by  the  wheel  speed  -  increasing  the  radial 
shift  at  the  tip  and  the  meridional  curvatures  -  reducing 
the  already  small  radial  shifts  near  the  hub. 
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The  increased  residence  time  and  radial  shift 
of  the  boundary  layer  particles  results  in  this  lou 
momentum  fluid  obtaining  a  larger  energy  transfer  from 
the  rotor  than  the  freestream  fluid.  This  nonuniform 
energy  transfer  results  in  circumferential  variations 
in  the  rotor  exit  total  temperature  field.  These  tem¬ 
perature  variations  are  largest  at  the  tip  of  the  blade 
due  to  the  influence  of  the  uheel  speed  on  the  radial 
shift  and  energy  transfer. 


6.  RECOMMENDATIONS  FOR  FUTURE  UQRK 


Much  of  this  project  was  devoted  to  the  devel¬ 
opment  of  the  calculation  technique.  Only  one  rotor  uas 
analyzed  in  detail  and  unfortunately  no  test  data  were 
available  for  this  rotor.  The  calculation  procedure 
should  be  applied  to  a  high  speed  rotor  for  which  reli¬ 
able  velocity  and  temperature  field  measurements  exist. 
There  is  a  problem  however  in  comparing  the  calculation 
results  with  measurements.  The  calculation  procedure 
stops  at  the  blade  trailing  edge  but  most  measurements 
are  made  some  distance  downstream  of  the  rotor.  The 
calculation  should  be  extended  to  include  some  reason¬ 
able  distance  downstream  of  the  rotor.  This  would  entail 
incorporating  a  boundary  layer  -  to  -  wake  transition 
model.  This  would  however  make  direct  comparisons  be¬ 
tween  calculated  results  and  measurement  data  possible. 

The  analysis  of  the  transonic  rotor  uas  re¬ 
stricted  to  velocity  ratios  created  within  the  blade  row. 
The  influence  of  inlet  casing  boundary  layers  or  up¬ 
stream  blade  wakes  entering  the  blade  row  was  not  con¬ 
sidered.  Such  an  analysis  should  be  made.  The  program 
as  it  exists  now  is  capable  of  making  such  an  analysis 
if  accompanied  by  a  pa r ti c le/f rees tream  mixing  or  inter¬ 


action  model 


Due  to  the  impact  the  boundary  layer  model  has 
on  the  resulting  trajectory  and  energy  transfer  calcu¬ 
lations  some  thought  should  be  given  to  incorporating  a 
more  elaborate  boundary  layer  calculation. 

Finally,  the  calculation  procedure  should  be 
extended  to  include  the  analysis  of  centrifugal  machines. 
The  radial  shift  considerations  are  a  little  different 
for  this  type  of  application .  Uhat  might  be  of  particu¬ 
lar  interest  is  modelling  the  inlet  casing  boundary 
layers  and  attempting  to  determine  if  there  is  a  rela¬ 
tionship  between  the  radial  migration  of  the  hub  casing 
boundary  layer  in  the  impeller  and  the  separation  location. 
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APPENDIX  I.  P  ROGRaI‘1  INPUT 


The  general  input  format  for  the  program  is 
given  in  Table  1.1  uith  sample  data  shoun  in  Table  1.2. 
The  program  operation  is  separated  for  absolute  velocity 
ratios  (lou  momentum  fluid  entering  the  blade  row)  and 
relative  velocity  ratios  (lou  momentum  fluid  generated 
within  the  blaae  row).  The  input  possibilities  for 
these  two  modes  will  be  explained  separately.  Common 
to  all  of  the  input  data  is  the  freestream  velocity 
field,  the  hub  the  shroua  contours,  ana  the  radial  lo¬ 
cation  of  the  particle  or  streamline  of  interest  at  the 
inlet  of  the  blade  rou.  The  velocity  ratio  and  flow 
angle  data  necessary  vary  uith  the  cases  considered 
(explained  below).  The  program  options  are  selected 
through  the  use  of  a  "trigger"  card  specifying  six 
triggers,  II,  12,  13,  14,  15,  and  16.  Trigger  13  de¬ 
termines  whether  absolute  velocity  ratios  (13=0)  or 
relative  velocity  ratios  ( 13  =  1 )  are  considered. 

I . 1  Absolute  Velocity  Ratio 

For  velocity  defects  in  the  inlet  flow  field 
it  is  necessary  to  specify  the  velocity  ratio  and  the 
relative  flow  angle  evolution  of  the  particle  through 
the  blade  row.  The  velocity  ratio  can  be  specified  as 
a  constant  or  as  some  function  of  the  axial  distance 
through  the  blade  row.  The  relative  flow  angle  evolution 


is  assumed  to  be  a  linear  variation  between  the  inlet 
and  exit  if  not  otherwise  specified  in  the  input.  Each 
of  these  cases  will  be  briefly  outlined  below. 

Case  1 .  ("trigger"  card:  0  □  0  0  0  0  ) 

For  this  case  the  velocity  ratio  is  an  input 
parameter  and  assumed  constant  through  the  blade  row. 

The  relative  flow  angle  evolution  is  assumed  linear 
from  the  inlet  angle  to  the  exit  angle. 

Case  2.  (  1  0  U  0  0  0  ) 

This  is  the  same  as  Case  1  except  now  the  velO' 
city  ratio  is  specified  as  a  function  of  the  axial  dis¬ 
tance.  The  relative  flow  angle  evolution  is  still  a 
linear  one. 


Case  3.  (  1  1  0  0  0  0  ) 

For  this  case  both  the  velocity  ratio  and  the 
relative  flow  angle  evolution  are  specified  as  functions 
of  the  axial  distance.  The  blade- to-blade  trajectory  of 
the  particle  can  now  be  arbitrarily  specified  allowing 
for  freestream/particle  mixing  or  interaction  models  to 
be  incorporated.  Note  that  the  inlet  relative  flow 
angle  is  obtained  from  the  freestream  flow  direction  and 
the  velocity  ratio.  Thus  the  relative  flow  angle  evolu- 


tion  should  be  specified  such  that  it  agrees  uith  this 
value  at  the  inlet. 


I . 2  Helative  Velocity  Ratio 

One  of  the  initial  goals  of  this  project  uas 
to  predict  the  radial  migration  of  blade  boundary  layer 
fluid.  This  mode  of  operation  can  consider  such  flous 
as  uell  as  freestream  (i.e.  out  of  the  boundary  layers) 
velocity  variations.  This  mode  of  operation  is  selected 
by  13  =  1.  For  relative  velocity  ratios  it  is  assumed 
that  the  relative  flow  angle  evolution  for  all  blade-to- 
blade  streamlines  is  (unless  otherwise  specified)  identi¬ 
cal  to  that  of  the  mean  streamline.  Thus  only  the  velo¬ 
city  ratio  needs  to  be  known  to  obtain  the  particle 
trajectory . 


Case  1 ■  (001000) 

The  relative  velocity  ratio  is  specified  in 
the  input  and  assumed  to  be  constant  through  the  blade 
row.  This  case  is  not  of  much  importance. 

Case  2.  (  1  0  1  0  0  0  ) 

Here  the  velocity  ratio  is  specified  as  a 
function  of  the  axial  distance.  The  relative  flow  angle 
evolution  is  that  of  the  mean  freestream  streamline. 


The  velocity  ratio  evolution  and  the  relative 
flou  angle  evolution  are  nou  specified  in  the  input  as 
functions  of  the  axial  distance  through  tne  blade  rou, 
thus  allowing  for  deviations  from  the  mean  streamline 
relative  flou  angle  evolution. 

Case  4.  (  0  0  1  Q  0  1  ) 

Up  to  this  time  it  has  been  assumed  that  the 
particle  follows  the  mean  streamline  through  the  blade 
row.  For  this  case  it  is  possible  to  specify  the  blade- 
to-blade  location  of  the  particle  of  interest.  This  is 
done  through  the  specification  of  Y/PITCH  in  the  input. 

A  value  of  0.0  corresponds  to  the  suction  surface,  0 .  “5 
to  the  mean  streamline,  and  1.0  to  the  pressure  surface 
of  the  blade  (still  assumed  out  of  the  boundary  layers). 
For  intermediate  values  of  Y/PITCH  an  interpolation  is 
done  between  the  surface  velocities  and  the  mean  stream¬ 
line  velocity  to  obtain  a  velocity  ratio  relative  to  the 
mean  streamline  velcity. 

Case  5.  (  0  0  1  1  0  0  ) 

For  this  case  it  is  possible  to  consider 
suction  surface  boundary  layer  particles.  The  desired 
streamline  in  the  boundary  layer  in  specified  as  a  y/5£ 


value  ranging  from  1 .00  for  the  streamline  coinciding 
uith  the  edge  of  the  boundary  layer  at  the  blade  exit 
(identical  to  the  freestream  surface  velocity)  to  0.0 
for  particles  being  "deeper"  uithin  the  boundary  layer 
(see  figure  6.).  The  velocity  defect  is  obtained  from 
the  boundary  layer  calculation  and  the  velocity  profile 
model  and  compared  to  the  mean  streamline  velocity. 

Case  5 .  (  0  0  1  0  1  0  ) 

This  is  iaentical  to  Case  5.  except  for  the 
pressure  surface  boundary  layer.  It  is  also  necessary 
to  input  the  number  of  points  involved  in  the  boundary 
layer  thickness  calculation  ana  the  fluid  dynamic 
viscosity  for  these  boundary  layer  cases. 


Table  1.1  General  Input  Format 


card  1 

card  2 

card  3 

card  4 

card  5 

card  6 

cara  7 

card  8 
< 

card  9 
card  10 
card  11 

ii 

it 

it 


(all  reals  read  in  F10.4  FORMAT) 
rotational  speed  (rad/s),  pitch  (m) 
inlet  temperature:  particle,  strm.  a,  strm.  B 
tip  LE  radius,  axial  location  (m) 
hub  LE  radius,  axial  location  (m) 
mid-merid.  distance  tip  radius,  axial  loc. 
mid-merid.  distance  hub  radius,  axial  loc. 
tip  TE  radius,  axial  location 
hub  TE  radius,  axial  location 
num.  strm Ins  -  IR,  num.  pts.  -  IZ  (214) 

-  begin  reading  freestream  - 

-  velocity  field  - 

1  -  IZ,  1  -  IR 

Read:  2  ,  R  ,L  ,  alpha  ,  beta  ,  uJss  ,  Ups 
from  Katsanis  Output 


then:  "trigger"  card  11,12,13,14,15,16  (612) 

:  Radial  location  of  entering  particle  on  interest 

then:  input  varies  depending  on  the  option  using 
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Sample  Data  Sets 


Set  1 

card  129: 
card  ISO: 
card  131: 
card  132: 
card  133: 
card  134: 
card  135: 
card  1 36 : 
card  137: 
card  138: 
card  139: 
card  140: 
card  141: 


"trigger"  card  (bI2) 

raaial  location,  %  span  (R) 

num.  pts.  in  velocity  ratio  distribution  (12) 
velocity  ratio,  axial  location 

"  ii  n  ii 

if  ii  it  it 

ii  ii  n  ii 

num.  pts.  in  flow  angle  distribution  (12) 
relative  flow  angle,  axial  distance 
ii  ii  ii  ii  ii 

'•  ii  ii  ii  ii 

ii  ii  n  ti  n 

start  next  data  set  uith  neu  span  location 


Set  2 


card  129: 
card  130: 
card  131: 
card  132: 


"trigger"  card  (612) 
radial  location,  ■/„  span  (14) 

Y/PITCH 

start  next  data  set  uith  neu  span  location 


Set  3 


card  129: 
card  130: 
card  131: 
card  132: 
card  133: 
card  134: 


"trigger"  card  (612) 

radial  location,  %  span  (14) 

num.  pts.  used  in  boundary  layer  calc.  (12) 

fluid  dynamic  viscosity  -  51  units 

y /£g  ,  b.  1.  streamline  of  interest 

start  next  data  set  uith  neu  span  location 


1  - 


APPENDIX  II.  PROGRAM  OUTPUT 


Examples  of  the  program  output  are  shown  in 
Table  II. 1.  The  headings  identify  the  mode  of  operation 
and  general  inlet  conditions.  The  particles  radial  lo¬ 
cation  and  velocity  is  calculated  at  12  axial  locations 
through  the  blade  row  at  increments  based  on  the  blade 
axial  chord  at  the  specified  inlet  radius.  The  axial 
and  radial  location  of  the  particle  is  given  at  each 
station,  as  well  as  the  velocity  ratio  and  the  velocity 
components  of  the  particle  and  the  freestream.  The  span 
(here  based  on  the  blade  height)  is  also  given  at  each 
axial  location  to  gain  a  quick  appreciation  for  the 
magnitude  of  radius  variations  calculated. 

The  blade  inlet  conditions  are  then  summarized 
for  three  different  streamlines,  the  particle  streamline, 
streamline  A,  and  streamline  B  (see  Figure  5.).  The 
velocity  components  and  flow  angles  are  given  for  each 
of  these  streamlines.  The  blade  exit  conditions  are  like¬ 
wise  given.  Lastly  the  energy  transfer  for  the  three 
different  streamlines  is  summarized.  Given  is  the  radial 
shift  of  each  streamline,  the  inlet  total  temperature, 
the  exit  total  temperature,  the  temperature  rise,  and 
the  exit  temperature  difference  between  the  particle 


streamline  and  streamline  B 
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Two  additional  parameters  are  included  in  the 
output  that  have  not  yet  been  explained.  They  are  labeled 
"Ij/VZB"  and  "E/P".  These  parameters  estimate  the  "cross- 
flow"  of  a  particle  and  the  resulting  displacement  in 
the  blade- to-blade  plane.  U/l/ZB  is  the  ratio  of  the 
component  of  the  particle  velocity  that  is  perpendicular 
to  the  relative  free stream  velocity  direction  and  the 
freestream  axial  velocity.  E/P  is  the  fraction  of  the 
blade  pitch  the  particle  will  be  displaced  as  it  flows 
through  the  blade  row.  Such  a  displacement  will  occur 
only  for  particles  that  have  at  some  time  a  relative 
flow  angle  different  then  that  of  the  freestream.  These 
parameters  were  included  mainly  out  of  the  curiousity 


of  the  author 
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